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ABSTRACT 

An ellipsoidal velocity distribution function is found as the result of maximizing the 
entropy with constraints on the number density, the mean velocity, and the stress energy 
tensor. The relation of this function to the macroscopic flow equations (continuity, 
momentum, etc.) is studied. The special case of one -dimensional channel flow Is ana* 
lyzed, It is shown that the ellipsoidal distribution function gives a consistent description 
of viscous How. That is, the parabolic flow profile and Newton's law of viscosity become 
natural consequences of the macroscopic flow equations , 



THE ELLIPSOIDAL VELOCITY DISTRIBUTION FUNCTION 
by N. Stankiewicz 
Lewis Research Center 

SUMMARY 

An ellipsoidal velocity distribution function is found as the result of maximizing the 
entropy with constraints on the number density, the mean velocity, and the stress energy 
tensor. The relation of this function to the macroscopic flow equations (continuity, 
momentum, etc.) is studied. The special case of one-dimensional channel flow is ana- 
lyzed. It is shown that the ellipsoidal distribution function gives a consistent description 
of viscous flow. That is, the parabolic flow profile and Newton's law of viscosity become 
natural consequences of the macroscopic flow equations. 


INTRODUCTION 

The concept of maximum entropy is well known in classical thermodynamics where 
it is used as a definition of the stable equilibrium condition. However, information 
theory, with its basic assumption that the disorder of a statistical system is a measure 
of its entropy, can provide a way of extending the maximum entropy concept to nonequi- 
librium situations. 

The ellipsoidal velocity distribution function is an example of a nonequilibrium (non- 
Maxwellian) distribution function that can be derived from a maximum entropy (disorder) 
principle (ref. 1) . Its name comes from the fact that the surfaces of constant density in 
velocity space are ellipsoids. 

The derivation of the distribution function rests on the supposition that the entropy is 
a maximum, subject to various constraints. The constraints are those moments of the 
distribution function which, it is hypothesized, are deduced from measurements made on 
the system. They are all local averages and can change in time and with position. An 
open (flowing) system is considered. 

The equations describing a flowing system are, in general, a hierarchy of equations 
involving the changes of a given moment of the distribution function in time; its transport 


into and out of an infinitesimal volume; the changes because of external forces; and 
finally the irreversible changes caused by particle interactions. 

This report is a study of the relation between the ellipsoidal distribution function and 
the macroscopic flow equations. The purpose herein is to show that the ellipsoidal dis- 
tribution function gives a consistent description of viscous flow. 

The procedure used in this report is in contrast to a previous work by this author 
(ref. 1) in which the relation between the Lagrange multipliers of different species within 
a single system (plasma) was investigated. 

Holway (ref. 2) and Kogan (ref. 3) show that the distribution functions obtained by the 
Grad moment method are linearized expansions of the more general functions obtained by 
using the maximum entropy procedure. That is, those terms which depart from being 
Maxwellian are linearized in references 2 and 3. This implies that the departure from a 
Maxwellian distribution function and, hence, from equilibrium is small. 

The present report takes the more general approach and leaves the distribution func- 
tion intact. The relation of the ellipsoidal distribution function to the macroscopic flow 
equations is therefore not restricted to near -equilibrium or to related assumptions con- 
cerning the ’'binariness*’ of collisions. 


ANALYSIS 

Maximum Entropy 

It is shown in reference 1 (and duplicated in appendix B of this report) that a dis- 
placed ellipsoidal velocity distribution function is the result of maximizing the entropy 
density 


S = 



f(ln f - l)d 3 v 


with the following moments of the distribution function as constraints: 
Number density: 


U) 


n ~ 



( 2 ) 


2 



Mean velocity density: 


n<v> 


■©y^ 


Stress -energy density tensor (per unit mass): 



fv v d 3 v 


(3) 


(4) 


(All symbols are defined in appendix A.) The nondimensional distribution function f is 
related to the density F of phase space by 


F = 



( 5 ) 


The moments defined by equations (2) to (4) are, by hypothesis (or by observation of 
the system), the independent variables of the system. That is, any other moment can be 
expressed in terms of these variables. 

Taking a variation of S (the variation operates only on the distribution function) and 
introducing the constraints by means of the Lagrange multipliers -O, -Jt, and -f, we 
obtain (as shown in appendix B) 

f = exp^-n-p • v - v ■ f • v) (6) 

where 

H = ~ 2( v ) ■ f (7) 

f" 1 * 2((vv) - <v><v>) (8) 


and 
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where 


A = 


det 


11 

r i2 

r !3 

12 

r 22 

r 23 

13 

r 23 

r 33 


( 10 ) 


The specification of the distribution function consistent with the assumed information 
is complete. That is, all the Lagrange multipliers are given in terms of the independent 
moments of equations (2) to (4). Equations (7), (8), and (9) define ju, f , and re- 
spectively, in terms of these moments. 


The Ellipsoidal Distribution Function and the Hierarchy of Flow Equations 

It is assumed in the preceding section that the constraints imposed by equations (2) 
to (4) are the results of measurements made on the system. However, these measured 
values may change in time, or from station to station within the system. There may also 
occur irreversible changes due to the internal interactions of particles. 

The changes in the constraints can be found by taking the appropriate moments of 
the Boltzman collision equation. Hence, corresponding to equations (2) to (4), for a sys- 
tem with no applied forces, the following equations can be derived: 


& + div p( v ) = 

at 

0 

(li) 

+ div p(vv) 

at 

= 0 

(12) 

Mh} + div p(v V V) - { 

at ' 

^P(vv}\ 

i at L 

(13) 


where p{=mn) is the mass density. 

The left hand sides of these equations are the changes due to time and to transport, 
while the right hand sides contain the irreversibilities (sources and sinks) of the flows. 
Equations (11) and (12) are, therefore, expressions of conservation of mass and momen- 
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turn. The trace of equation (13) is the energy equation whose right hand side would in- 
clude wall interactions. No use is made of this equation in this report. 

An important point to notice about this set of equations is the occurrence of the next 
higher flow average - namely, p(vvv) . Also, the equation for p(vvv ) would contain 
p(vvvv) . Seemingly a complete solution to any flow problem thus entails the solution of 
an infinite set of equations for all the flow averages. The usual resolution of this diffi- 
culty is to truncate arbitrarily the hierarchy of equations. The point of truncation usually 
coincides with the exhaustion of names for the higher terms - the heat flux generally 
being the last term considered. 

The appearance of higher flow averages is not a problem if the distribution function 
is known, because the needed term is then calculable. 

The maximum entropy technique outlined in the Maximum Entropy section (see p, 2) 
provides the most probable distribution function, consistent with the constraints con- 
sidered. All higher moments are then calculable and expressable as functions of the in- 
dependent constraining moments. This procedure is essentially an a priori truncation 
method for the set of equations (11), (12), and (13). From appendix B we get 

P<Vj v k> = P[< v i>< v j v k> + < v i v j>< v k> + < v jX v i v k> _ 2 < v i> < v j> <v] (14) 

Equation (14) is used in the next section to simplify equation (13) for analyzing a 
special case of viscous flow. 


The Ellipsoidal Distribution Function and Viscous Flow 


The purpose here is to show that the ellipsoidal distribution function gives a correct 
description of physical reality, in that it contains the elements of viscous flow. 

The case considered is steady-stale incompressible channel flow in which the flow 
velocity is one«dimensional in the x direction, that is, 


<v> = (v ) e 


X / X 


and the pressure tensor is 


P = p[<vv> - <v)<v>] = 


P wy 0 


P_ r P 0 

xy 


(15) 


(16) 


0 


0 


p 
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where P is the scalar pressure. 

Equations (11) and (15) with the conditions of steady-state and incompressible flow 
immediately specify that 


(v x > = V(y) 


( 17 ) 


Because of equations (15) and (11), equation (12) can be rewritten as 


div p((vv) - <v><v>) = 0 
Then, with the use of equation (16), equation (18) becomes 


(18) 


3P. 


0x 9y 


23U o 


(19a) 


ap 


= 0 


9y dx 


(19b) 


It is assumed that the scalar pressure is constant in the y- directi on, that is, transverse 
to the flow. Equation (19) then becomes 


dP. 

+ 

dx dy 


dp ■ ~~EL b o 


( 20 ) 


where P = P(x) and P xy = P xy (y) . 

The final flow equation to be considered is equation (13), which expresses the 
changes in the stress -energy tensor but contains the higher moment p(v i v-Vj c ). Nor- 
mally, this would be an unknown term, and for the analysis to continue, one would need 
to make some assumption concerning it. However, because of the maximum entropy 
procedure used in this report the most probable distribution function is known and this 
moment is calculable and given by equation (14). 

Equation (13), with the use of (14), then reduces to 


P 


dV 


dy 



( 21 ) 
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If the interaction term is linearized by 



( 22 ) 


it is seen that equation (21) contains Newton's law of viscosity - namely. 



(23) 


where 77 = P fv is the coefficient of viscosity, and v is an interaction frequency, or the 
reciprocal of a characteristic relaxation time. 

Equations (20) and (23) can now be solved for the velocity profile. Each term in 
equation (20) must be equal to a constant so that for the case of constant viscosity equa- 
tion (20) becomes 


-K 

dx 


(24a) 


d 2 V _ K _ r 
dy 2 v 


(24b) 


where equation (23) was used in equation (24b). The scalar pressure, according to equa- 
tion (24a), then varies linearly with the distance down the channel. And, solving equa- 
tions (24), with the conditions 


V = V for ^ = 0, y = 0 
° dy 


V = 0 for y = ±a 


J 


(25) 


yields the parabolic velocity profile 


a 2 


( 26 ) 
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It is thus seen that the ellipsoidal distribution function leads to a meaningful physical 
result. 


CONCLUDING REMARKS 

The principle of maximum entropy (disorder) consistent with the moments n, n{v), 
and n(v v) as constraints leads to the ellipsoidal velocity distribution function. The con- 
straints are introduced into a variational equati on by means of Lagrange multipliers and 
these unknown multipliers are subsequently related to the constraints. The constraining 
moments are a set of independent macroscopic variables of the system. 

In general, the flowing (open) system is described by a hierarchy of flow equations 
(continuity, momentum, etc.) which shows how the moments can vary with respect to po- 
sition or time. However, the flow equations are not a closed set in that they always in- 
volve an unknown higher moment. 

Because the constrained moments constitute an independent set of variables, the un- 
known moment is expressible in terms of the constraints. That is, all moments are cal- 
culable because the distribution function is known. Moreover, the distribution function 
is the most probable one consistent with the constraints considered. Thus, the principle 
of maximum disorder provides a consistent truncation for a given set of flow equations. 

For the specific case of the ellipsoidal distribution function and for one-dimensional 
incompressible channel flow it is shown that the truncation of the flow equations contains 
Newton's law of viscosity. And this leads to a parabolic velocity profile. 

An extension of the method to include heat flux and other higher moments is desir- 
able and further studies might yield criteria predicting the point when higher moments 
must be introduced to describe adequately the flowing system in general. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, November 12, 1968, 

129-02-08-05-22. 
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APPENDIX A 


SYMBOLS 


A 

c 

F 

f 

h 

A 

I 

J 

K 

m 

n 


( W T z> 
(x, y, z) 

fi 

a 

r 

v 

M 

r 

r 

p 

n 

O) 


determinant of the F matrix, defined in eq. (10) 
constant 

density in velocity -position (6 -dimensional) space 
velocity distribution function 
Planck's constant 
unit tensor 

Jacobian, determinant of transformation matrix 

constant 

mass 

number density 

scalar pressure; or, with subscripts, component of pressure tensor 
pressure tensor 
entropy 
mean velocity 

velocity vector; components 
velocity vector defined in eq. (B15) 
components of position vector 
parameter defined in eq. (B16) 
matrix of Lagrange multipliers 
coefficient of viscosity 

vector whose components are Lagrange multipliers 

interaction frequency 

velocity vector defined in eq. (B12) 

mass density 

Lagrange multiplier 

degeneracy factor 
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Subscripts: 

i,j,k running indices (= 1,2,3) 

Mathematical symbols: 
det determinant of 

div divergence in coordinate space 

grad gradient in coordinate space 

grad y gradient in velocity space 

-1 inverse of 

dot or scalar product 
* vector 

5 tensor symbol 

A average value of 
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APPENDIX B 


IDENTIFICATION OF THE LAGRANGE MULTIPLIERS 

The variational equation used in deriving the ellipsoidal distribution function is 

SS - ftSn - pt • 6^n(v}^ - f . . 6^n(vv)j = 0 (Bl) 

Substituting equations (1) to (4) into equation (Bl) yields 

-go — J 5f(ln f+Sl+p-v+v - f • v)d 3 v = 0 (B2) 


Because the variation §f is arbitrary, the integrand must vanish so that 

( -Jl-ju- v -v • f • v 
f = e 


(B3) 


which is the ellipsoidal distribution function in equation (6) . 

There are many ways to arrive at equation (7) and (8), but the simplest is as follows. 
Taking the gradient in velocity space of f gives 

gra<y = - f(p + 2v • f ) (B4) 


in which the symmetry of F is used in equation (B4). 

When integration is performed over all velocities, the left hand side of equation (B4) 
vanishes because f must vanish for infinite velocities: 


gradyf d 3 v = ^ ^ = € i f [f(°°) " f(-°°)]d 2 v = 0 


Hence, equation (B4) becomes 


0 = 



2<v> • 


r) 


(B5) 
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where equations (2) and (3) were used. Finally, 


M = -2{v) • r 


(B6) 


which is equation (7) of the main text. 

Substituting equation (B6) into equation (B4) and forming the diadie give 

(grad v f)v = 2ff * (<v)v - vv) (B7) 

and noticing that 

grad^fv) = if + (grad v f)v (B8) 

one arrives at 

grad v <fv) - If = 2fr • (<v)v - vv) (B9) 

If equation (B9) is integrated over all velocities, the gradient term will again vanish. 
After simplifying one gets 


!= 2F • ({vv) - <v><v>) (BIO) 

Equation (BIO) immediately defines equation (8); that is, 

2({vv) - <v><v>) (Bll) 

The integration of equation (2) to get the explicit form of n in equation (9) proceeds 
as follows. First, transform to the relative velocity variable which is given by 

v = (v) + f (BX2) 

From this it follows that 

d 3 v = d 3 £ (B13) 

Then, with the use of equation (B6), the number density of particles in an incremen- 

tal volume of velocity space d 3 v becomes 
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The easiest way to carry out the indicated integrals is to complete the square of 
each component {| j, ! 2 , i; 3 ). The following nonorthogonal transformation will accomplish 
this result: 


where 



1/2 

r !2 

r i3 


W 1 


1 11 

r l/2 
1 11 

1/2 

A 11 

*1 

W n 


0 

/« \l/2 

/ P 33 \ 

^33 

^2 

2 


Vii) 

2 r ll%3 ^ 

W 3 


0 

0 

( a y /z 

\ 33/ 

h 


(B15) 


~\ 


^33 “ r ii r 22 " r i2 


^23 “ 2 ^ r i2 r 13 “ r il r 23^ 


a det r - r 11 r 22 r 33 + 2 r 12 r 13 r 23 _ r n r| 3 - r 22 r 13 - r 33 r 2 2 


(B16) 


It follows from equation (B15) that 


,3 

d w = Jd ij 


= A V2 


(B17) 


where J is the determinant of the transforming matrix in equation (B15). The integra- 
tion in equation (B14) can now be carried out: 



n . • r- <v> 

_ w ^ 3 jr 3/2 A -l/2 e -S2+(v> ■ f- <v> 


/ 


2,3 

-w d w 


(B18) 


which is equation (9) . 

Because of the independence of the Lagrange multipliers, one is able to compute any 
integral of f through a process of differentiation of n with respect to the Lagrange 
multipliers. 

Therefore, note that 

n = j f d^v = n(n, ju, f) (B19) 


and 


3n 
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fv i d°v = - n{v.) 


(B20) 


Likewise, since 


then 


<V 4 ) . i /iv, 


^T‘ = '^/ fV >= C| 3 v / fV i d 3 v 'n/ fV i V k " <V < T i> - < v iV <B21 > 


k n 


Lastly, 


!v i v j > = ; / fv i v j d3v 


14 


i 



and 


a(v,v,) 

= < v i v j> <V “ < v i v j v k> < B22 ) 

k 


Equations (B22), (B21), and (B20) can now be rewritten in terms of first, second, 
and third order differentiations of n, K this is done, then equations (B20) and (B2 1) 
verify the results given by equations (B6) and (Bll). 

Equation (B22) is the defining equation for (v^.v^) . Substituting equation (Bll) 
yields 


< T i v J v k> = < v i v j> <V 


®< v i v j> 

94 k 


= <Vj>< v k> 



< v i><V) 


= <vj> <V 


< Vi > 


3<Vj) 

9M k 



9<Vj) 

84 k 


(B23) 


Now using equation (B21) gives 

<Vj\> = < v i v j> <V k > + <Vk>< v j> + < v jV ^ v i> ' 2<v i ><v j )<v k ) (B24) 

Multiplying equation (B24) by p then gives equation (14). 
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